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Abstract — Optimal  radar  waveform  design  for  target  detection 
has  been  addressed  in  prior  research  literature  under  various 
assumptions  regarding  noise  and  clutter.  A  common  model  of 
the  radar  scene  in  work  of  this  kind  is  a  linear  time-invariant 
(LTI)  operator  with  additive  Gaussian  noise  that  acts  on  the 
transmitted  signal  to  produce  the  received  signal.  This  model 
is  intrinsically  ill-suited  to  dynamic  scenes  or  moving  radar 
platforms  because  it  cannot  account  for  Doppler.  This  paper 
introduces  scene  models  based  on  Hilbert-Schmidt  class  (HS) 
operators  on  the  space  of  finite-energy  signals.  This  category  of 
models  generalizes  the  LTI  category  in  the  sense  that  every  LTI 
operator  is  also  a  HS  operator,  but  the  HS  class  includes  opera¬ 
tors  that  account  for  frequency  shifts  as  well  as  time  shifts  and 
are  thus  suitable  for  modeling  radar  scenes  involving  Doppler. 
Every  HS  operator  is  uniquely  expressible  as  a  superposition  of 
elementary  time  and  frequency  shift  operators,  thus  providing  a 
convenient  interpretation  of  a  scene  in  terms  of  these  physically 
meaningful  operations  on  the  transmitted  signal.  Application  of 
this  perspective  to  waveform  design  for  target  detection  in  noise 
and  to  optimal  receiver  processing  for  a  given  waveform  for  target 
detection  in  clutter  and  noise  are  demonstrated. 

I.  Introduction 

The  idea  of  regarding  a  radar  scene  as  a  linear  or  nonlinear 
operator  that  transforms  transmitted  waveforms  to  received 
signals  has  been  used  in  several  radar  contexts.  In  particular, 
this  perspective  has  been  applied  to  optimal  waveform  design 
for  detection  in  noise  and  clutter  under  various  resource 
constraints  (e.g.,  [5],  [6]).  Often,  the  scene  is  modeled  as  a 
linear  time-invariant  (LTI)  system.  This  class  of  systems  is 
very  well  understood  in  signal  processing  and  is  of  value  as 
a  radar  scene  model  for  static  scenarios,  but  it  is  inadequate 
for  situations  involving  Doppler.  Specifically,  since  sinusoids 
are  eigenfunctions  of  LTI  systems,  this  mathematical  model 
offers  no  means  for  the  scene  to  impart  a  Doppler  shift  to  a 
narrowband  transmitted  waveform. 

The  primary  goals  of  this  paper  are  to  examine  a  broader 
class  of  linear  operators,  the  Hilbert-Schmidt  class,  in  the  role 
of  radar  scene  models  and  to  illustrate  how  they  can  be  used 
in  the  design  of  waveforms  and  optimal  receiver  processing 
for  target  detection  in  noise  and  clutter. 

II.  Hilbert-Schmidt  Operators 

Denote  by  L2(R)  the  Hilbert  space  of  complex  signals  /  : 
ffi.  — >  C  with  inner  product 

if,  9)=  f  f(t)g*{t)dt 

JR 


and  norm 

11/11  =  (</,/>)1/2  = 

A  LTI  system  is  a  linear  operator  S  :  L2(M)  — >  L2(M)  that 
can  be  expressed  as  convolution  with  a  fixed  h  £  L2(M);  i.e., 

Sf(t)  =  [f*  h](t)  =  f  f{r)h(t  -  r)  dr 

J  R 

Under  fairly  general  conditions,  limiting  arguments  permit 
extension  of  the  domain  of  S  to  allow  its  application  to  ideal 
(Dirac)  impulses  f(t)  =  d(t  —  to)  and  complex  sinusoids 
/(f)  =  eif2t,  even  though  these  are  not  elements  in  L2(R). 
The  function  h  is  known  as  the  impulse  response  of  the  system 
because  h  =  S6.  If  /  is  a  complex  sinusoid  with  frequency 
O, 

Sf(t)  =  J  eiQjTh(t  -  t)  dr 
=  eiQt  f  h{u)e~mu  du 

JR 

=  hmm 

The  response  of  the  system  S'  to  a  sinusoidal  input  at  a  given 
frequency  is  thus  a  sinusoid  at  the  same  frequency  whose 
amplitude  and  phase  are  determined  by  the  complex  constant 

h(Q)  =  (  h(u)e~lQu  du 

Jr 

(i.e.,  the  Fourier  transform  of  h  evaluated  at  frequency  Q), 
which  is  known  as  the  frequency  response  of  the  system  at  this 
frequency.  The  point  in  reiterating  these  well  known  properties 
of  LTI  systems  here  is  to  emphasize  that  such  systems  cannot 
manifest  a  shift  in  the  frequencies  present  in  their  input  signals. 
They  are  thus  inadequate  for  modeling  dynamic  radar  scenes 
involving  Doppler. 

The  Hilbert-Schmidt  (H-S)  class  [2]  is  a  larger  collection 
of  linear  operators  on  L2(R)  that  can  accommodate  Doppler 
shifts.  Consider  the  displacement  operator  D(a,  r)  :  L2(R)  — > 
L2(M.)  defined  by 

D{a,r)f{t)  =  elat  f(t  —  r) 

i.e.,  the  operator  D(a,  r)  shifts  the  signal  f  by  t  in  time 
and  by  a  in  frequency.  A  Hilbert-Schmidt  class  operator  S  : 
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L2(R)  — >  L2(K)  is  one  that  can  be  written  as  a  superposition 
of  displacement  operators  [4] 

S  =  /  /  s(a,r)D(a,T)  dadr  (1) 

Jr  Jr 


where  the  “coefficients”  s(a,r)  form  a  square-summable 
function  of  the  variables  a  and  r;  i.e.. 


|s(ct,r)|2  da  dr  <oo 


The  more  usual  characterization  of  the  H-S  class  is  as  the 
linear  operators  on  i2(K)  satisfying  Tr(5Sb)  <  oo.  In  terms 
of  the  expansion  fl) 

Tr(S'S^)  =  f  f  \s(a,r)\2  dadr  <  oo 

Jr  Jr 

The  relation  (1)  can  be  seen  as  an  invertible  transformation 
between  L2(R2)  (i.e.,  the  space  of  square-summable  functions 
of  the  variables  a  and  r)  and  the  H-S  class.  The  inverse 
transform  is 

s(a,r )  =  Tr  (SD(a,r)^) 


Note  that  restricting  (1)  to  a  =  0  and  writing  h(r)  =  s(0,  r) 
gives 

Sf(t)=  f  h(r)f(t  —  t)  dr  =  [f  *  h](t) 

Jr 


So  LTI  operators  are  a  subclass  of  the  H-S  class. 


A.  Scene  Modeling 

Consider  the  monostatic  radar  situation  in  which  a  target 
consists  of  a  discrete  collection  of  point  scatterers  with  un¬ 
known  position  and  unknown  Doppler.  Invoking  narrowband 
and  far-held  assumptions  [7],  the  H-S  operator  description  of 
such  a  scene  is 

S=  Y.  st(a,T)D(a,T)  (2) 

(a,r) 

where  the  sum  is  taken  over  all  pairs  of  Doppler  shift  a  and 
time  delay  r  corresponding  to  a  constituent  point  scatterer 
in  the  target  and  the  complex  weight  st(a,r)  represents  the 
strength  of  the  point  scatterer  at  delay-Doppler  position  (a,  r). 
If  the  scatterers  comprising  the  target  are  clustered  around  a 
nominal  “target  center”  (an,  To)  in  time  delay  and  Doppler,  it 
is  possible  to  re-write  (2)  in  the  form 


S  = 


22  st(a,T)D(a,T) 
(a,r)G  iV(0,0) 


D(a  0,T0) 


In  this  expression,  N (0, 0)  is  a  neighborhood  of  (0, 0)  in  (a,  r) 
coordinates  measured  from  the  target  center.  The  contributions 
of  the  scatters  within  this  neighborhood  are  summed,  each 
with  a  displacement  operator  D(a,  r)  capturing  its  differential 
range  and  Doppler  relative  to  the  target  center.  The  operator 
D(a o,t0)  shifts  the  response  of  this  entire  cluster  of  point 
scatterers  by  ao  in  Doppler  and  by  To  in  time,  corresponding 
to  the  range  and  radial  velocity  of  the  target  center  relative 


to  the  radar.  This  model  extends  to  multiple  targets  by  direct 
superposition. 

Similarly,  clutter  arising  from  collections  of  point  scatterers 
may  be  expressed  in  the  form 

c  =  22  sc(a,T)D(a,r) 

(a,r) 

The  nature  of  such  clutter  may  be  modeled  through  assump¬ 
tions  on  the  weight  function  sc  :  R2  — »  C;  e.g.,  deterministic, 
Gaussian  or  Markov  random  field,  or  other  statistical  char¬ 
acterizations.  In  the  stochastic  case  there  are  advantages  to 
assuming  a  continuum  model  for  the  clutter 

C=  /  /  sc(a,  r)D(a,  r)  da  dr 

Jr  Jr 

where  sc  G  L2(R2). 

Suppose  that  the  clutter  is  characterized  as  a  random  field 
on  R2 .  That  is,  the  function  sc  is  considered  a  random  field 
which  is  taken  to  have  zero  mean  and  covariance  function 
R((a,r),  (a',T')).  As  such,  sc  has  a  Karhunen-Loeve  expan¬ 
sion  [1,  p.  71]  in  terms  of  the  eigenfunctions  of  the  covariance 
operator 

OO 

sc(a,T )  =  "^2  \  X,Ti°Jn-T] 

3= 1 


where  E(£j)  =  0  and  E (£j£k)  =  $jk  for  j,  k  =  1,  •  •  •  ,  oo.  Here 
the  (j>j  are  the  orthonormal  eigenfunctions  of  the  covariance 
operator;  i.e., 


R((a,T),(a' ,T'))<f>j(a' ,t')  da' dr'  =  A j<f>j(a,T) 


and 


/  /  <j>j(a,  r)(j>k(a,  r)  da  dr  =  Sjk 
Jr  Jr 

Defining  the  clutter  eigen-mode  operators  1’,  by 

Tj=  /  (j>j(a,T)D(a,T)  dadr,  (3) 

Jr  Jr 


they  satisfy 

Tr^rl)  =  6jk 


and  the  clutter  operator  can  be  written  as 

OO 

C  =  E  \/Wf  w 

3= 1 


In  practice,  only  a  certain  number  of  the  eigenvalues  will  be 
appreciable  and  the  sum  can  be  truncated.  Note  that  this  clutter 
model  does  not  assume  any  form  of  stationarity  for  the  clutter. 

Stochastic  target  models  can  be  developed  in  a  similar  way 
in  terms  of  a  localized  random  field.  With  such  a  model,  S 
would  take  the  form  of  a  stochastic  linear  combination  of 
orthonormal  target  eigen-mode  operators  A  j  with  associated 
eigenvalues  jij 


Sk 


M 

'22 

3  =  1 


D(a  0,T0) 


where  M  is  an  appropriate  truncation  point. 
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III.  Detection  in  AWGN 

Consider  a  monostatic  radar  with  baseband  waveform  w  G 
L2(R).  At  the  carrier  frequency  Cl,  the  waveform  is  D(Cl,  0)ru. 
Employing  a  H-S  scene  model  with  one  illuminated  point 
scatterer  at  range  r  and  radial  velocity  v  relative  to  the  radar, 
the  return  is  s(ao,  To)D(ceo,  ro)[-D(fI,  0)u>]  where  tq  =  2 r/c 
and  ao  ~  2 vCl/c  (c  is  the  speed  of  light).  The  process  of 
mixing  the  received  signal  back  to  baseband  yields 

x(t)  =  D(—Cl,  0)s(ao,  r0)D(a0,  t0)D(CI,  0 )w(t) 

=  e~znT°eiaots(a0,To)w(t  -  r0) 

=  s(ao,T0)e~inTo  D(a0,T0)w(t) 

For  an  extended  target  in  zero-mean  additive  white  Gaussian 
noise  n{t),  (2)  gives 

x(t)  =  ^  e~lQTst{a,T)D(a,T)w(t)  +n(t) 

(a.r)  (5) 

=  Tw(t)  +  n(t) 

where  the  H-S  class  target  response  operator  T  :  L2(R)  — > 
L2(K)  is  defined  by 

T  =  ^2  e~lQTst(a,T)D(a,T) 

(«>t) 

To  test  the  hypothesis  Hi  that  the  received  signal  x  has 
the  form  (5)  against  the  null  hypothesis  Hq  that  x(t)  =  n(t), 
observe  that  under  Hi,  x(t)  is  Gaussian  with  mean  Tw(t)  and 
variance  er2  (i.e.,  the  variance  of  the  noise  process  n).  Under 
Hq,  x(t)  is  Gaussian  with  mean  zero  and  variance  a2.  Due  to 
the  Gaussian  nature  of  this  problem,  detection  performance 
(e.g.,  in  the  Neyman-Pearson  sense)  [3]  is  optimized  by 
maximizing  the  signal-to-noise  ratio,  which  is  equivalent  to 
maximizing 

||Tw||2  =  ( Tw,Tw )  =  (T^Tw,  w) 

where  T'  denotes  the  adjoint  of  T.  Hence  detection  perfor¬ 
mance  is  optimized  for  a  given  signal  energy  by  choosing 
the  waveform  w  to  be  an  eigenfunction  of  the  non-negative 
definite  Hermitian  operator  P  =  T'T  corresponding  to  its 
largest  eigenvalue. 

If  the  constraint  is  power  rather  than  energy  (i.e.,  with 
|w(f)|  <  B  for  all  t),  it  is  generally  possible  to  improve  on 
this  eigenfunction  solution.  This  analysis  will  be  presented  in 
a  later  paper. 

IV.  Detection  in  clutter  and  noise 

As  in  the  preceding  section,  assume  a  target  has  response 
operator  T  which  is  in  the  H-S  class  in  L2(M)  and  that  the 
noise  is  white  and  Gaussian  with  variance  cr2.  In  addition,  the 
transmitted  signal  w  is  affected  by  a  H-S  class  clutter  operator 
C.  The  signal  presented  to  the  receiver  is  thus 


“decision  instant”  in  [5]).  This  value  is  given  by  r(to)  = 
( HTw  +  HCw  +  Hn)(to)  and  has  associated  signal  to 
interference  plus  noise  ratio 


SINR(fo) 


\HTw(t0)\2 

E\HCw(to)  +  Hn(to)\ 2 


(6) 


As  in  [5],  the  initial  objective  is  taken  to  be  maximizing  the 
SINR  at  time  to,  given  w,  by  suitable  choice  of  H.  Ultimately, 
once  the  optimal  H  is  known  as  a  function  of  w,  the  waveform 
that  maximizes  SINR  when  paired  with  its  corresponding 
optimal  receiver  processing  H  will  be  sought.  This  second 
optimization  step  is  only  discussed  superficially  in  this  paper, 
however. 

As  a  H-S  operator  on  L2(R),  H  may  be  represented  as 
an  integral  operator  with  kernel  <I>  :  M2  — ►  C;  i.e.,  for  any 

/  e  £2(R), 

Hf(t)  =  I  dr 

Jr 


Consequently, 

E\Hn(t)\2 


=  E 


<£(£,  r)n(r)dr J  <£(£,  u)n(u)du 

=  a2  J  \$(t,r)\2  dr 


Denoting  hto  =  <l>(Vo,  ■),  the  properties  of  the  transform  kernel 
imply  hto  G  L2(M),  and  the  preceding  expression  gives 


E|tfn(to)|2  =  a2||/Uo||2  (7) 


For  the  clutter  term, 

E\HCw(to)  |2  =  J  J  ht0(r)ht0{u)Cw(r)Cw{u)  dr  du 

=  ( ht0,Gcht0 )  (8) 

where  Gq  is  a  non-negative  definite  Hermitian  operator  de¬ 
fined  in  terms  of  the  waveform  w  and  the  clutter  operator  C 
by 

Gcf(t )  =  /  E[Cw(t)Cw(r)\  /(r)  dr 
Jr 

The  presence  of  the  expectation  operator  in  this  expression 
indicates  that  the  clutter  operator  may  be  regarded  as  random. 
If  this  is  the  case,  the  ensuing  analysis  is  unaffected  provided 
that  the  noise  n  has  zero  mean  (as  already  assumed)  and  is 
independent  of  C.  For  deterministic  (known)  clutter 

Gc  =  CPWC f 


s(t)  =  Tw(t)  +  Cw(t)  +  n(t) 

The  receiver  processes  s  with  a  H-S  operator  H  and  makes 
its  detection  decision  on  the  basis  of  the  statistic  Hs(to),  the 
response  of  this  operator  at  a  particular  instant  to  (called  the 


where  Pw  is  the  rank-one  projection  operator  from  I2(W)  onto 
the  one -dimensional  subspace  spanned  by  w;  i.e.,  for  /  G 

L2(R)  . 

Pwf  =  (f,w)  W 
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In  the  stochastic  case  the  expansion  (4)  gives 


Gc  —  E  57  \/ A? Xk^jG; r.j  Pw r  1 


u'=i  fc=i 


OO 

5>,i  Vvrj. 

i=i 


Finally,  the  signal  term  in  the  numerator  of  (6)  may  be 
written  in  inner  product  form  as 


\HTw{t0)\2  =  {ht0,PTwht0)  (9) 


where  Ptw  is  the  rank-one  projection  operator  from  L2(R) 
onto  the  one-dimensional  subspace  spanned  by  Tw.  Substitut¬ 
ing  (7),  (8),  and  (9)  into  (6)  yields 


SINR(fo) 


(ht0 1  PTwht0) 

(htn,(a2I+Gc)ht0) 


(10) 


where  I  denotes  the  identity  operator  on  L2(R).  In  this 
form,  maximization  of  SINR  may  be  treated  as  a  generalized 
eigenvalue  problem.  The  operator  Pq  in  the  numerator  is 
positive  definite  and  of  rank  one.  Since  Gc  is  non-negative 
definite,  the  operator  ( cr2l+Gc )  in  the  denominator  is  positive 
definite  as  long  as  a 2  >  0.  Thus  SINR  is  maximized  by  setting 


ht0  —  fcmax  —  (cr~I  +  Gc)  XGTw 


To  calculate  the  maximal  SINR  value,  denote  G  =  (a2 1  + 
Gc)  and  observe  that  G  is  positive  definite  and  self-adjoint. 
Thus  /imax  =  G~1!‘1Tw  and  the  maximal  value  of  SINR(fo) 
is 


SINR 


max 


(G~1GTw,  PTwG~1/2Tw) 
(G-WTw,  G1/2Tw) 


Denoting  f  =  Tw  and  using  the  self-adjointness  of  G  gives 


SINRmax 


(G-1/2/,  PtwG-1!2  f) 
(G-1/ 2/,Gi/2/) 

(/>  G~1GpTwG~1G  f) 

(/./) 


The  maximal  value  of  this  ratio  is  exactly  the  maximal  eigen¬ 
value  of  the  non-negative  definite  operator  G~xGpTwQ~xl2 _ 
Since  this  is  a  rank-one  operator,  its  trace  is  equal  to  its  one 
positive  eigenvalue.  Hence, 


SINRmax  =  Tr  (G-^Pt^G-1/2) 

=  (G-^jjG-1/2/) 

=  (f,G-1f)  =  (Tw,G~1Tw) 
=  (Tw,  (ct2I  +  Gc)~1Tw) 


Note  that  this  agrees  with  [6,  Eq.  (11)]  when  the  noise  is  white 
and  for  the  special  case  of  LTI  operators. 

The  waveform  design  problem  is  reduced  to  finding 

■u’max  =  argmax  (w,  T\a2I  +  Gc)1Tw)  (11) 

wGW 

where  IE  is  a  permissible  class  of  transmit  waveforms.  Note 
that  this  is  a  non-linear  problem  because  the  operator  Gc 


depends  on  w.  When  the  clutter  operator  G  is  known  (de¬ 
terministic)  then 


(aH  +  Gc)-1  =  I- 


GP,„Ct 


{w,(CW  +  a2I)w) 


and  so  the  design  problem  becomes 


Wmax  =  arg  max 
w€W 


^(w,T*Tw) 


|  (w,  T^Cw)  |2  \ 

(w,  (GtG  +  a2l)w) ) 


V.  Conclusions  and  Future  Work 


This  paper  has  proposed  modeling  of  radar  scenes  by 
Hilbert-Schmidt  class  operators  on  C2(K),  arguing  that  this 
class  of  models  extends  the  commonly  used  LTI  class  of  op¬ 
erator  models  to  accommodate  both  time  delay  and  frequency 
shifts.  Some  remarks  on  stochastic  clutter  modeling  within  this 
framework  were  also  presented. 

If  a  target  with  known  operator  model  is  to  be  detected  in 
additive  white  Gaussian  noise,  the  optimal  waveform  subject 
to  a  total  energy  constraint  has  been  observed  to  arise  as  an 
eigenfunction  of  an  operator  formed  from  the  target  model 
operator.  In  the  presence  of  a  known  clutter  operator,  the 
receiver  processing  that  maximizes  the  SINR  at  a  chosen  time 
instant  was  shown  to  be  characterized  by  the  solution  to  a 
generalized  eigenvalue  problem. 

Future  work  will  investigate  optimization  of  the  waveform 
in  the  presence  of  clutter  and  subject  to  an  energy  constraint. 
Investigation  of  the  use  of  this  operator-theoretic  perspective 
to  optimize  joint  waveform  and  processing  design  under  power 
rather  than  energy  constraints,  as  mentioned  in  Section  III,  is 
also  underway. 
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